\ 

N76-21063 

RADIO DETECTIONS OF COMETARY MOLECULAR TRANSITIONS: A REVIEW 
L. E. Snyder 


I . li'TTRCDUCrKU 

The various radio searches for comet ary rnolecviles finally came to 
fruition with the radio detections of CH, CH, HCN and QigCN and several 
unidentified species in Comet Kbhoutek (1973f) and the detection of H 2 O 
emission from Conet Bradfield (1974b). As a result of these detections we 
have learned scxnething about wiiy past radio searches of other ccxnets for 
molecules were less than successful. Also, it is now possible to identify 
and discuss some of the similarities and differences between conetaiy mole- 
cules and interstellar molecules, particularly with regard to the excitation 
and chemistry. Finally, given the observed projected densities and resulting 
gas production rates, the feasibility of future radio n»lecular observations 
of conets can be discussed. 

II. THE IDENTIFIED RADIO MOLECULES 
All of the identified radio molecules have been found in two ccxrets, 
Kohoutek (1973f) and Bradfield (1974b), and are listed in column 1 of taole 
1. Colvram 2 gives the quantum numbers of each detected transition and the 
corre^onding rest frequency is in column 3. Column 4 lists the antenna 
tenperature for each line and column 5 shows whether the radio line was 
detected in emission (E) or absoiption (A). The antenna half -power beam- 
width (HPBW) is in column 6 and the reported date of the observing period is 
in column 7. Column 8 has the approximate ratio of the heliocentric 
distance R to the geocentric distance A for the period of observation. 

Column 9 gives the calculated projected density N (molecular density per cm~2 
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Table 1 

Radio Detections of Identified Cometary Molecules 


Molecule 

Transition 

v(MHz) 

Ta(K) 

a/e 

HPBW 

Date 

R(AU)/A(AU) 

N(can~^) 

Reference 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 





OaVET 

KOHDOTEK 

(1973f) 

/ 



CH3CN 

J,K=6,0-5,0(Vg) 

110,712.2 

0.6 

E 

67"x63" 

1973, 1&5 Dec 0.9/1. 4 

3x10^®* 

Ulich cind 





NRAO 


' 


Conklin (1974) 


6,3-5,3(Vg) 

110,709.6 

0.4 

E 

36-ft. 





HCN 

J=1-0,F=2-1 

88,631.8 

0.8 

E 

75" 

1973,15-16 

Dec 0.5/1. 2 

1.5x10^3* 

Huebner , 






NE«0 




Snyder & Buhl 






36-ft. 

1974, 3-6 

Jan 0.3/0. 9 


(1974) 

CH 

J=1/2,F=1-1 

3,335.5 

0.2 

E 

14" 

1974, 2-7 

Jan 0.3/0. 9 

4x10^“*** 

Black et al . 






fCO 

25.6-m 




.(1974) 

OH 

J=3/2,F=2-2 

1667.4 

0.09+ 

A 

3.5'xl9* 

1973 Dec 

<0.7/l,2> 

3.3x10^^ 

Biraud et al. 


P=l-1 

1665.4 

0.07+ 

A 

^Jancay 



17 

4,1x10'^-^ 

(1974) ~ 


F=2-2 

1667.4 

0.09+ 

E 

Meridian 

1974 Jan 

<0.75/0.8> 



P=l-1 

1665.4 

0.05+ 

E 










CXMBT 

BRADFIEU) 

(1974b) 




H^O 

■^K-K+=®16~^23 

22,235.1 

0.15 

E 

1.5' 

1974,9 May 

1.22/1.15 


Jackson, Clark 





I'lERDC 

120-ft 




and Donn(1975) 



* total projected density calculated assumirg Boltzmann statistics 

** projected density for the ground state lambda doublet 

+ peak antenna tenperature 



integrated along the line-of-sight) and the detection references are listed 
in colvtrnn 10 , 

Several interesting points are shewn by table 1. Three of the conetary 
molecules were detected in the millimeter- wavelength region (CHgCN at 2.7 ran, 
HCN at 3.4 ran and H 2 O at 13.5 ran) while the remaining -two were detected in 
the older, more developed centimeter wavelength region (CH at 9 cm and CH 
at 18 can) of the radio spectrum. Thus the millimeter wavelength region is 
proving to be very usefiiL for the detection and study of cemetary molecioles. 
Given the noise limitations of currently available millimeter and centimeter 
radio receivers, all of the detected cemetary radio lines had antenna taiper- 
atures (see <x)limm 4) which were close to the detection limit and were 
weak when conpared with mast of the cinrent interstellar moleciiLar observa- 
tions. This suggests that for future comet searches succjessful molecular 

12 

detections will require a projected density in excess of about 10 molecules 

n 

per cm"*^ averaged over the half-power bearavidth. Of cx>virse if the antenna 
beam is larger than the cometary molecular cloud, the projected density re- 
quired for detection in the optically thin case will incjrease in proportion to 
the square of the ratio of the beam diameter to the molecular cloud diameter. 
Thus the typically smaller beamwidths in the millimeter spectral region (see 
cioluran 6 ) cxjnsiderably reduce the beam filling factor and allow the detection 
of polyatomic cemetary molecules which have low projected densities. Success- 
ful radio molecular observations of Kohoutek were made for R less than 1 AU 
(see colvnnn 8 ), when the molecular production rate was high, and for A no 
greater than about 1.4 AU so beam dilution was hot excessive for the molecu- 
lar species of interest. '* 

The detections of CH 3 CN (Ulich and Conklin 1974), HCN (Huebner, Snyder 
and Buhl 1974) and H 2 O (Jackson, Clark and Donn 1975) were the first radio 
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observations of expected parent or mother molecules in ccmets (Potter and 
Del Duca 1974; Huebner 1970) but because of the possibility of nonequilib- 
rium excitation, the projected densities and resulting production rates for 
these molecules must be interpreted with caution. For exanple, CH 3 C 3 Si was 
detected and observed with good signal-to-noise only in the vibrationally 
excited Vg = 1 state vdiich is approximately 50QK above the zero^xjint eneigy. 

It should be noted that vibrationally excited CHgCN never has been detected 
in the galactic molecular clouds; the only rraiotely ccnparable mLcrowave 
transiticais observed in galactic sources are the emission lines from rotation- 
al transitions of the vibrationally excited SiO maser (Snyder and Buhl 1974; 
Davis et 1974; Thaddevis et 1974; Buhl et al. 1974) ^\hich originate in 
circumstellar envelopes of late-type stars, i^proximately two weeks after the 
013 ( 2 ^ Vg = 1, J = 6-5 transitions were observed by Ulich and Conklin (1974), 
Buhl, Huebner and Snyder (1975) searched lor both Vg = 1, J = 5-4 and ground 
vibrational state J = 5-4 lines of CHgQJ. The detected CHgCN lines were 
strong enough to confirm the Ulich and Conklin (1974) detection but the 
signal-to-noise ratio was so low that the new data were not siifficient to 
serve as an independent identification of canetary CHgCN (see Figure 2 in 
Buhl et 1975). The ground state lines did not have the intensities which 
would be ejqjected from a Bolt 2 mann distribution. Hence both the hi^ exci- 
tation tenperatiare required by Vg = 1 and the transient nature of the signal 
strength suggest that the CH3CN Vg = 1 excitation was time-varying. There- 
fore the projected density (column 9, table 1) and production rate of 10^® 
obtained by assuming a Boltzmann distribution over vibrational states of 
CH3CN probably are too large and a more meaningful production rate is 
'v 10^7 _ i()28 g-1 obtained frc*n the weak ground state signal (Buhl ^ 1975) . 

Evidence of time-varying emission also was found for HCN. The data of Huebner 
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et al. (1974) show multiple HOST Doppler cxinponents . A caiponent with zero 
shift showed increases as well as decreases in the daily average while other 
conponents appeared and disappeared during the observations. Because there 
was no observed evidence for nonequilibrium excitation of HGn and because 
the observations lacked spatial resolution, the calculated production rates 
of 1.2 X 10^^ s"^ (for r^ = 10“^ km) to 3 x 10^^ s~^ (for r^ = 3 x 10^ km) 
are based on two assunptions: a Boltznann distribution and isotropic 

outstreaming of the gas (Snyder, Huebner and Buiil.1975). Ihe difficulties 
of calculating an H 2 O production rate from the radio observations of Canet 
Bradfield (Jackson et 1975) are due to both the high excitation (447 aii”^) 
required to populate the 6^0 level and the relatively long radiative lifetine 
for the 6 ig - 533 transition. As a result, HgO production rates fomd from 
applying Bolt 2 raann statistics to radio obseiwations tend to give numbers 
that are too high and a more realistic production rate probably is around 
10^^ (see table 1 of Jackson et (1975)). 

Radio observations of cometary CH prcxnise to be useful for studying 
radiative puirping models of CM and hence may be directly applicable to 
current research on galactic CM source excitation. In 1973 Decanber, 

Biraud ^ (1974) detected CM absorption in Qxnet Kohoutek at 1665 and 

1667 MHz. Turner (1974) used the NRAO 140 ft. telescope to confiim their 
observations while the comet was still approaching perihelion. After 
perihelion, however, Biraud et al . (1974) found that the CM reappeared in 
onission and tliat ultraviolet punping of the ground state lambda doublet 
could roughly account for the observed CM intensities as a fvuiction of time. 
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III. THE UNIDEINTiriED RADIO MOLECULES 
Two groips of radio observers foimd imidentified lines wiiile searching 
Comet Kohoutek for other molecules. Giguere and Clark (1975) were iising the 
NRAO 140-f t . telescope to search for the 8190 MHz line of CH when a weak 
imidentified emission feature at 8188.82 MHz was detected. J.K.G. Watson 
(1974) has suggested that this line may be one ccnponent of the 6 ]l 6 “ ^23 
rotational transition of the NH 2 radical. Buhl, Huebner and Snyder (1975) 
were using the NRAO 36-ft. telescope to search for vibrationally excited 
SiO when fairly strong unidentified lines at 86,247.1 and 89,010.5 MHz were 
found. They present arguments that the 86,247 line can not be assigned 
to nearby transitions of ethanol, acetone or vibrationally excited silicon 
monoxide and suggest that one or both lines may belong to unstable species 
which are decay products of more cotiplex molecules. 

IV. NEGATIVE RESULTS 

' In recent years numerous conetary searches have been conducted at 
several radio observatories but not all of the negative results have been 
formally reported. The reported negative results (which I am aware of at 
present) are ordered by frequency in table 2. Column 1 lists each molecule 
and indicates when a vibrationally excited state was sought; column 2 gives 
rotational qiiantum numbers and column 3 the corresponding rest frequency. 
Column 4 lists the ipper limit reached for the peak-to-peak antenna toiper- 
ature T^ except when noted otherwise; column 5, the radio telescope used; 
coliinn 6 the iialf -power beamwidth (HFBW) in arc minutes; and column 7, the 
reference . 

It shovild be noted frcxn columns 4 and 5 that very low antenna toiper- 
ature limits were reached for H 2 O, 


H^OO, NHg and CI^(CN )2 using favorable 
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Table 2 

Reported Negative Results 


Molecule 

Transition 

V (MHz) 

T^(K) 

Telescope 

HPBW 

Reference 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 




(30MET 

BEUNETT (1969i) 



H 2 CX) 

^ir^io 

4,829.7 

0.30 

NRAD 140-ft 

6.6 

Huebner and Snyder (1970) 

H 2 O 

^16"523 

22,235.1 

2.5 

NRL 85-ft 

2.3 

Clark ^ (1971) 




COMET KOHOOTEK (1973f) 



H 2 CO 

^ll"ll 0 

4,829.7 

0.03* 

MPI 100-m 

2.6 

Schroder et al. (1974) 

OH 

l/2,F=l-0 

4,765.6 

0.03* 

II II 

" 

II II 

OH 

5/2,F=3-3 

8,189.6 

0.20 

MEAD 140-ft. 

3.6 

Giguere and Clark (1975) 

HC 3 N 

1-0,F=2-1 

9,098.3 

0.15 

II It 

3.2 

II II It 

HC 3 N( 2 \) 7 ) 

1-0,F=2-1 

9,156.1 

0.25 

II M 

II 

II It II 

HCN(V 2 ) 

6 , ^doub. 

9,423.3 

0.40 

II It 

3.1 

II II II 

H 2 O 

®16"^23 

22,235.1 

1.07 

MPI 100-m 

0.7 

Churchwell et al. (1975) 


II 

II 

0.02 

APD 26-m 

1.4 

Avery and Andrew (1974) 


II 

II 

nr** 

NEBOC 120-ft 

1.5 

Jackson et al. (1975) 


3,2 

22,834.2 

0.03 

ABD 46-m 

1.4 

Avery and Andrew (1974) 

CH 2 (CN )2 

^16“'^07 

23,084.2 

0.04 

AKD 46-m 

II 

II II II 

NH 3 

1,1 

23,694.5 

0.7 

MPI 100-m 

0.7 

Churchwell et al. (1975) 

CH 3 OH 

42 ~ 4 i 

24,933.5 

0.75 

II . II 

II 

II II 

CH 3 C 2 H 


85,457.2 

0.4 

NRAD 36-ft. 

1.3 

Buhl ^ al. (1975) 

((313)20 

^ 20"^11 

86,222.9 

0.4 

II II 

II 

II tl 

Sio(V=l) 

2-1 

86,243.3 

0.1 

If II 

II 

II II 

HNCO 

‘*04“^03 

87,925.2 

0.3 

M II 

II 

II II 

HNCO 

'*13-3i2 

88,239.0 

0.4 

II II 

II 

II tl 

HCN( 2 v 2 ) 

1-0,F=2-1 

89,087.9 

0.3 

II II 

II 

II tl 

X-ogen(HCO'^) 1-0 

89,188.6 

0.3 

II II 

II 

II II 

"HNC" 

1-0 

90,665 

0.3 

II It 

It 

II tl 

HC 3 N 

10-9 

90,979 

0.4 

II II 

II 

II II 


1-0 

113,491.0 

0.2 

HEAD 36-ft. 

1.1 

Ulich and Conklin (1974) 

CO 

1-0 

115,271.2 

0.5 

II II 

II 

II II II 


* brightness tanperature limits 

** not reported 
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beaniwidths without achieving detections in Ccmet Kohoutek. The well known 
anall error in the topocentric ephemerides used by most radio observers has 
been discussed by several authors (e.g. Avery and Andrew 1974) and probably 
it had little effect on the outcome of most observations. However, ephoneris 
accuracy is most essential for ccnpact soiirces and the high excitation re- 
quired by the - 623 H 2 O line suggests that in most ccaiets it will be found 
only in a rather oonpact region. Thvis, both the weak intensity of the H 2 O 
emission reported by Jackson ^ (1975) fron Ccmet Bradfield and the 

generally low elevation of Ccmet Kohoutek suggest that atmospheric H 2 O atten- 
uation coi 5 )led with a small pointing error may have been responsible for the 
negative H 2 O results found for Comet Kohoutek. 

The upper limit for H 2 CO emission found by Schroder st (1974) for 
Ccmet Kohoutek was at least 100 times more sensitive than the i?>per limit 
found by Huebner and Snyder (1970) for Comet Bennett because of the narrower 
beam width used for the Ccxret Kohoutek search. Detection of the 111 - iio 
H 2 CO transition in ccmets may not be as straightforward as it appears, how- 
ever, because if the excitation conditions are similar to the typical inter- 
stellar case, then the H 2 CO punping mechanism may force the absorption mode 
to predominate. The CH observations and radiative punping calculations of 
Biraud ^ (1974) suggest that preperihelion searches for H 2 CO absorption 

may be more successful than postperihelion searches in future comets, 

V. SUMMARY 

We have learned that successful radio spectroscopy of ccmets almost 
always reqiiires going to the detection limits of current instnjmentation , 
which corresponds to a projected nunt>er density (in the initial energy level 
of the transition) of approximately 10^^ molecules per cm^. Many of the 
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raxiio observations appear to be consistent with a straightforward appli- 
cation of the fluid d 5 mamic model (Huebner and Snyder 1970; Clark ^ 

1971) with the possible exception of H 2 O (as discvissed by Jackson et 
(1975)) and other moleciiles reqxiiring unusual excitation for raxiio detection. 
Thus for most molecules of interest, reasonable estimates for the success of 
future searches can be obtained by starting fran an expected production rate 
estimate, possibly based on the Ccmet Kohoutek results, and working backward 
to see if a projected density above 10 ^^ cm~2 can be reached without 
excessive beam dilution. 

Several of the successfxil radio detections have denonstrated that not 
all radio lines lead to a physically meaningful production rate due to the 
excitation conxlitions required for ctetection. Hence, as in the study of 
galactic molecular soxnrces, non-LTE mechanisms will haye to be introduced 
to explain the radio observations of both high lying transitions (e.g. H 2 O 
and vibrationally excited CH 3 CN) and of easily inverted transitions (e.g. 

CH and possibly H 2 CD) . As a result , we can now start thinking about using 
future ccxnetary observations to qpiquely test several of the molecular 
punping models which have been proposed for galactic molecular sources. 

The unidentified lines found in Ccmet Kohoutek were at least as strong 
as the identified lines nearby in the spectrum and have not , as yet , been 
detected in the galactic molecular sources. This suggests that the radical 
apd molecular ion chonistry of comets should be given more attention in the 
radio spectrum. In fact, it may be more interesting if radio astroncmers in 
the future perfomed a few frequency scans on comets instead of limiting 
thonselves to grinding away on known transitions of stable molecules. 

Finally, as a practical consideration, we have learned that both a good 
set of radio ephenerides and (coordination among various groxps of raxiio 
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observers are essential to success. It is doubtful that I would have as 
many conetary radio detections to discuss in this review had it not been 
for the efforts of Drs. T. A. Clark (NASA), D, K. Yeomans (Con^), Sci. Coip.), 
B. G. Marsden (SAG) and others in preparing and distributing useful ephener- 
ides especially designed for radio observers. Ihe exceptional coordination 
between various radio observers and observatories optimized the use of 
valuable radio telescope time and was chie mostly to the work of Drs. W, E. 
Howard, (I®AO) , S. P. Maran (NASA), and R. W. Hobbs (NASA). Vv'ithout their 
efforts, many of the radio molecular observations of Ccxnet Kohoutek would 
not have occurred. 

I wish to thank the many radio observers \dio sent me their comet 
results for use in this review and acknowledge partial support for this work 
from NSF grant OP-34200 to the University of Virginia. 


241 



REFERENCES 


Avery, L, W., and Andrew, B. H. 1974, Astr.J. , 79, 1322. 

Biraud, F. , Bourjois, G. , Crovisier, J., Fillit, R. , Gerard, E. , and Kazes, I. 
1974, Astr. and Ap. , 34, 163. 

Blaxds, J. H. , Qiaisson, E. J., Ball, J. A., Penfield, H., and Lilley, A. E. 

1974, Ap. J. (Letters) . 191, IA5. 

Buhl, D., Huebner, W. F. , and Snyder, L. E. 1975, this pxiblication. 

Buhl, D. , Snyder, L. E., Lovas, F. J., and Johnson, D. R. 1974, J. (Letters), 
192, L97. 

Churchwell, E., Landecker, T. , Winnewlsser, G. , Hills, R. , and Rahe, J. 1975 
These proceedings, p. 297. 

Clark, T. A., Bonn, B., Jackson, W. M., Sullivan, W. T. Ill, and Vandenberg, 

N. 1971, Astr. J. . 76, 614. 

Davis, J. H., Blair, G. N., Van Till, H., and Thaddeus, P. 1974, Ap. J. (Letters), 
190, U17. 

Giguere, P. T. and Clark, F. 0. 1975, Ap. J. (in press for 15 June). 

Huebner, W. F. 1970, Los Alamos Report IA-4542-MS, 

Huebner, W. F. , and Snyder, L. E. 1970, Astr. J. , 75 , 759. 

Huebner, W. F., Snyder L. E., and Buhl, D. 1974, Icarus , 23, 580. 

Jackson, W. M. , Clark, T. , and Bonn, B. 1975, this piJblication. 

Potter, A. E., and DelDuca, B. 1964, Icarus, 3, 103. 

Schroder, R. , Wendker, H. J., and Stunpff, P., 1974, Astr . and Ap . 37, 417. 

Snyder, L. E., and Biahl, D. 1974, Ap. J. (Letters) , 189, L31. 

Snyder, L. E., Huebner, W. F., and Buhl, D. 1975, Comet Kohoutek Workshop , 
NASA/MSFC , (NASA Special Publication, in press). 

Thaddeus, P., Mather, J. , Davis, J. H. , and Blair, G. N. 1974, J. 

(Letters) , 192, L33. 

Turner, B. E. 1974, Ap. J. (Letters) , 189, L137. 

Ulich, B. L., and Conklin, E. K. 1974, Nature , 248, 121. 

Watson, J. K. G. 1974 (preprint). 


242 



DISCUSSION 


H. Keller : I would like to ask if you would define a flare, and also 
do you have any idea what the scale length of H^O was ? 

L. Snyder : Okay. Your second question, I think should be directed to 
Jackson, who’s sitting here. 

The first question - I can define a flare as anything that excites a vibra- 
tionally excited state, hi this case we think it may have been a knot that started 
passing from the head of the comet to the tail about the 1st of December. But 
you can take the usual definition of flare and see what it does in the radio spec- 
trum, or you can take the radio definition and then say if this excitation is un- 
usual, it’s certainly flaring in the radio. 

L. Biermann; Could you say anything concerning the upper limit for the pro- 
duction rate of CO given by Conklin and Ulich? It’s dipole moment is very 
small. 

L. Snyder ; They gave upper limits for CN and CO. First of all, I’d just 
like to say a word about the upper limits. Given the current state of telescope 
technology and receiver technology, if you can’t see a molecule, the upper 
limit on the projected density for the state of interest usually comes out to be 
about 10^^ or 10*^ , This in turn, given reasonable cometary dimensions and 
a ground-state rotational transition, will invariably lead to a production rate 
on the order of 10 or 10 . But what I think happened in the CN upper limit, 

was that they got a beautiful spectrum of vibrationally excited CH 3 CN, two K 
components. And if I were to guess, I would say that if one could remove the 
methylcyanide lines, one might see a velocity broadened CN line. 

L. Biermann; What about CO ? 

L. Snyder: A production rate less than 10^^ would be consistent with 
the CO negative results. 

Voice : Is that a ground state transition ? 

L. Snyder : Yes. 1-0. 

F. L. Whipple : I want to congratulate Dr. Snyder and all the workers here, 
but I wanted to remove one misapprehension, at least in my opinion. I think 
that’s off the record that we consider cometary astronomy in a very bad state. 

(Laughter. ) 
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DISCUSSION (Continued) 


L. Snyder : You mean it's come full circle. 

F. L. Whipple ; It hasn't gotten very far yet, 

(Laughter. ) 

W. Jackson; I wanted to make a few comments about the radio observations. 
One is that the comet Bradfield production rates, as we realized, are completely 
ridiculous, in terms of the water production rate. I think that it points up the 
fact that what we may be seeing is a reflection of the excitation mechanism. 

Even the production rate for methylcyanide as calculated originally by Ulich and 
Conklin, and I think that's the original number, is too low if you use reasonable 
values of the radius of the comet - I mean I think they use something like 10 
kilometers or 20 kilometers for the radius of comet Kohoutek, and I think some- 
thing more reasonable is like 5. The signal goes as the square of the radius of 
the nucleus, so even that suggests that indeed methylcyanide must be produced 
either in a flare or in a nonequilibrium mechanism, as the vibrational excitation 
indicates. 

The beam dilution, it turns out, is really a function of the lifetime of the 
cometary species; the beam dilution in formaldehyde is greater than the beam 
dilution in water, or it should be. 

L. Snyder; It's also a function of wavelength and telescope diameter, so if 
you put the proper configuration - 

W. Jackson; Yeah. Okay, okay, I agree. — But I'm just saying the over- 
all cloud size is going to be dependent upon the photodissociation rate. And it 
turns out for water the photodissociation rate ought to be almost a factor of 4 
lower than the photodissociation rate of formaldehyde. 

H. Keller ; What is the photodissociation rate for H 2 O? 

W. Jackson ; The lifetime I get is something like 2 times 10'^ seconds. And 
for formaldehyde it is 5 or 6 times 10^ seconds. 

H. Keller ; At one AU? 

W. Jackson : One AU. Everjrthing I say is referred to lAU. 

I've puzzled over why we saw water in Bradfield rather than water in 
Kohoutek. If I may go to the board for a minute, I think part of it may be due to 
the way we were pointing in Bradfield. 
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DISCUSSION (Continued) 


I don't know much — I'm not really a radio astronomer, but I was up there 
for both of the observations, and in Kohoutek, if this is the earth — 

We were pointing the telescope at low elevation through a large atmospheric 
mass. The comet was down close to the horizon all the time. 

Bradfield was much higher - near the pole, and was circling around the pole. 
Now, the noise in the case of Kohoutek increases measurably compared to the 
noise for Bradfield. 

The effectiveness of the integration is much better for Bradfield because 
even though we may have integrated for 4 or 5 hours on Kohoutek most of it 
was under conditions where I would say that the signal to noise ratio was poor- 
er than under the 4 or 5 hours for Kohoutek. 

L. Snyder; Okay. If you combine this, a r secant Zfor the atmosphere cor- 
rection, say, with even a quarter of a beam width mispointing, the signal could 
disappear in the noise very easily. I'd just like to add a little bit to that, and I 
think that the nonequilibrium detections — methyl cyanide, water vapor, HCNj- 
I think are going to prove to be very important for understanding what happens in 
the interstellar clouds. Because if you consider for a minute — if you get enough 
telescope time and enough observational time, as the French did on the OH, as 
you track the comet, and as you can see signal variations with heliocentric dis- 
tance, it gives you very direct measures, then, on what the UV pumping, for 
example, is doing to a particular molecule. You just can't do it in the inter- 
stellar medium. 

I think it's a very promising area of future research, which has come 
out of this. 

G. Herzberg: May I ask to what extent you guard — you have to guard against 
accidental overlapping of some background source. I mean your comet, while 
you are averaging for several hours, goes over a lot of the sky. Now, there may 
by chance be — 

L. Snyder: That's right. That's right. 

/ 

Now, the primary way we guard against this is by looking at each scan as 
it comes in, before we do any stacking, do very careful editing. The secondary 
way we guard against it is to make sure that we are clear of the known sources, 
and then a third way, if any detection is suspect, the common way is just to go 
back and point in that direction after the comet's gone by. And if there's a source 
there, then the signal will come in even though the comet's gone. 
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E. Ney ; Maybe this is crazy, but what makes infrared so easy is that 
you can find a comet in the daytime, peak up on the signal, get the telescope 
tracking, and then go integrate on hard wavelengths. 

What are the chances that you’ll find something bright in the radio to track 
on? For example, broad band thermal emission from the dust? 

L. Snyder ; There were several thermal emission experiments tried 
at various frequencies, and I think the only one that was successful was the inter- 
ferometer experiment of Hobbs and Maran and Webster and — who else was on 
it — at Greenbank. And when I talked to Bob Hobbs a few minutes ago, they were 
still squabbling over the results, because again, it’s a very, very weak signal. 
It’s less than 70 flux units, Bob says. 

E. Ney ; It’s just not bright. 

L. Snyder ; There’s nothing bright, given present receivers. Now, if we 
can cut receiver noise maybe in half, maybe the time will come when we could 
do that. But given the beam dilution and given the receivers, it's a pretty im- 
jxassible job. We can’t even map now, given our beams and our signal to noise. 

We’d like to reach the place where we could map on comets, too, which is 
a sensitivity problem. 

So as far as I know — Dave, did you have something to add to that? 

D. Buhl: There was — I believe it was the French group detected the comet 
at 1 millimeter in the continuum. If the receivers are improved considerably, 

1 millimeter might be the frequency for that type of thing. 

P. Wehinger: In regard to the background noise that you were commenting 
on, in the case of Bradfield, it was going through the galactic plane as a back- 
groimd noise. Are there any problems with this compared to the case where 
Kohoutek was well out of the galactic plane ? 

L. Snyder; Why don’t you ask the man who owns one, sitting in front of 
you? Jackson did the Bradfield observations. 

W. Jackson; Well, the way the radio observations are done, we are 
looking first where we expect the comet to be 10 minutes later. And then 
we’re looking where the comet is now. So if it was in a galactic noise 
source, the background signal measured ten minutes earlier would have 
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included the galactic source also. Subtracting the background would 
eliminate the galactic source. 

L. Snyder: All right. You’re beam switching there? 

W. Jackson : Yes. 

B. G. Marsden : Pd just like to distinguish between errors in the ephemeris 
and errors in the use of the ephemeris. 

(Laughter. ) 

How accurately do you actually want the ephemeris ? 

L. Snyder : Well, okay. I brought along some numbers. It seems to me that 
if we're going to probe the nucleus of the comet, we want to probe at millimeter 
wavelengths . 

This means that our largest beam width will be on the order of 80 seconds 
of arc, and our smallest one, for all practical purposes, given today's tele- 
scopes, will be about 40 seconds of arc. So we don't want anything probably 
more accurate than 10 seconds of arc, because that's the pointing accuracy 
limit, but if you coidd take it down to the pointing accuracy — 

B. G. Marsden : Down to 30 seconds, say. 

L. Snyder: Well, 10 seconds would be great — 


B. G. Marsden: Then down to 10 seconds. 

L. Snyder: Right. 

B. G. Marsden: The custom is, of course, that epheme rides are published 
by combining the geocentric coordinates of the earth for universal midnight, or 
ephemeris midnight, with the heliocentric coordinates of the comet for that time. 
Now, this is absolutely useless, because nobody ever observes that, but this is 
the convention. What you really want, I suppose, is for us to antedate the helio- 
centric coordinates of the comet by the light-time, and produce an astrometric 
ephemeris. 

But even then you'd be getting into trouble with parallax, if you aren't 
already when you're getting down to it at ten seconds. But this is the convention 
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of the lAU; I suppose we have to ask Commission 20 if it will change its con- 
vention. 

L. Snyder ; Yes, well, times are changing in radio, so 1 think it would be 
a fair request, because as Jackson points out, when you have a low-lying comet, 
and you're trying to do water, and in addition you have a half a beam width error, 
as this pair of slides showed, it can just take the signal right down into the noise. 
And in fact, there was some hint of a signal on water in Kohoutek; but it really 
couldn’t be proven, couldn’t be verified above the noise. 

Isn't that correct? 

W. Jackson: Yes, There are some measurements that hint that there might 
have been some water present. But not enough that you could look at it and say 
with confidence, especially since I'm not a radio astronomer. 

L. Snyder: And this optical checking business probably isn’t any better than 
a beam width where it's tracked with the telescope mounted beside the radio 
telescope, because the axes don’t exactly coincide. 

B. G. Marsden: Well, we have provided in a few cases on request, special- 
ized ephemerides for certain people, actually for the location of their radio 
telescope. How useful is it to do this after the event? 

L. Snyder; Then it’s usually too late. Except for the cases where we have 
a strong signal of an imidentified feature and we want to get the rest frequency 
on it and determine whether it's ethyl alcohol or acetone or what have you. 

B. G. Marsden ; Did you have ethyl alcohol, perhaps? 

L. Snyder: I have an aimouncement from the Washington Post here that I 
could read, but I won’t. 

(Laughter. ) 

W. Jackson ; I’d like to make one more comment about the radio measure- 
ments in the nonequilibrium excitation mechanism. If you examine the infrared 
measurements of Meisel as I originally saw them — and I'm quoting from a 
preprint — the original CN to OH concentration was such — by the infrared 
measurements — you would predict that there was more CN in the comet than 
OH. 
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And this depended, as I understood the analysis, on a thermal mechanism 
for the production of the CN infrared emission. That contradicts the UV meas- 
urements as we know them, of CN and OH. We understand the excitation mech- 
anism of the CN and OH emission in the UV region. 

All I’m trying to say, really, is that one of the exciting things about the 
possibilities of observing comets in the infrared and radio regions, it may lead 
to more information about the chemistry of comets, since you are measuring 
transistions that may be excited by new mechanisms. 

M. Mumma r I'd like to say a few words on that, which bear also on the 
question of the accuracy of the ephemerides. This winter a group of people 
here at Goddard, and also at NRAO, including Dave Buhl, myself, Ted Kostuk 
Steve Cohen and Peter von Thuna, from Arthur D. Little, built a new instrument 
called an infrared heterodyne spectrometer, which promises the ability to 
measure IE emissions from parent molecules in the inner coma, with 5 arc 
second resolution in the neighborhood of 10 microns. 

This is obviously going to be an important new tool for study cometary 
processes, and it’s one that is now available. 

So I would say to Brian that we could use ephemerides that reflect the 
pointing accuracy that we can now achieve with this instrument. 

Voice ; I hate to sound skeptical, but I just wonder what criterion you use 
for the identifications you have in the microwave spectra. For example, can 
you tell me how many random noise spectra you would feed in for each identi- 
fied spectrum? In other words, is the chance 1 in 10, 1 in 20 of this being 
spurious ? 

L. Snyder : In our own particular data, I think that the vmidentified spectra 
stand alone, in that we have significantly no scans in which they’re not showing up 
above the noise leveL So I don’t think we have a statistical problem there. 

In the HCN data that we report, our common practice is to break the scans 
into groups, stack the groups and see if we can still see the line, and then try 
mixing the groups from different days, and see if we can still see the line. In 
the case of the central features in the HCN, which are all above the noise, we 
can. 

Now, in the case of the methylcyanide data (both the Vg = 1 state and the 
ground state data) if the features that were appearing at about 2 to 1 or 2. 5 to 1 
did not match the velocity which had been reported, say, a couple weeks earlier 
with about 4 to 1 signal to noise by Ulich and Conklin, then we couldn’t even 
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report them. I show them here as a confirmation of their detection, but they 
do not qualify as independent of their detections. 

But i don't really think that's the point in identifying a radio line. Or if it 
is, it doesn't explain enough. In order to identify a radio line, first of all, as 
we've discussed earlier, you have to make sure the pointing is on the comet 
and not on some galactic object. And then you have to go through a series of 
local oscillator shifts, and see if the line moves as you shift the local oscillator. 
If it moves in your filters or your autocorrelator, because you're moving the 
box around an absolute spectrum, so as you shift one way the line should move 
one way, and as you shift the other way the line should move the other. Any real 
feature that's coming in, say, from the comet or a galactic source, will move, 
and any noise feature will not move. And when you build up signal to noise on a 
feature that moves several times. The statistical argument, then, boils down 
to what the proper intensity is in terms of signal noise. 

E. Ney; I'm afraid this may even be almost insulting, but I assume that you 
handled differential refraction. You have to have your ephemeris, of course, 
but you have to also put in refraction. 

L. Snyder : Well, we take that out with our pointing. 

That's integral to our pointing procedure. 

E. Ney: But differential refraction, as you track down toward the horizon — 
what do you do about that ? 

L. Snyder: Well, as we track down toward the horizon, we try to correct 
our pointing accordingly. 

E. Ney : Toward the differential refraction? 

L. Snyder : We know where we're trying to point, and we know 
where our objects are that we have as standard check objects. I think that ques- 
tion has come up before on Bennett and I think that's the best answer I can give 
you for now. 

E. Ney: But it's big. You go another air mass and it's another minute of 
arc, so if you're pointing at an ephemeris position you'll be out of the beam. 

D. Buhl: (inaudible) 


250 



DISCUSSION (Continued) 


L. Snyder: What Dave Buhl said was that it's in the telescope pointing pro- 
gram, and it's a problem that you have with any kind of source. It's essen- 
tially taken care of in a series of pointing sections every few months and it's 
programmed in. 

Does that answer your question ? Okay. 

Because — well, yes, that's a good point. We track galactic sources over 
the horizon, for example, qxiite commonly, and I must say that the horizon 
pointing, which takes into account differential refraction, sometimes is much 
better than the zenith pointing. We sometimes get better signal to noise as the 
source is going over the horizon than we do at zenith on the 36-foot telescope. 

So I'm happy that the corrections are proper from the data that I've seen 
in the past. 

W. F. Huebner : Two models are being proposed for the formation of comets. 
One of them deals with the formation in the solar nebula. The other one, in a 
companion nebula, the composition of the comet would therefore be different in 
both of them. 

If the formation of the comet were to occur in the solar neb\ila, you would 
expect chemistry which is more like equilibrium. I put the emphasis on "more 
like," not on "equilibrium." 

Can you tell us what the ratios are of the "exotic molecules" to H 2 O in inter- 
stellar space? Because I think that would be an important criterion to differen- 
tiate between the two models. 

L. Snyder : Okay, I can give you some numbers. But one thing to keep in 
mind, which I pointed out before, that given the models, probably the only mole- 
cule that we have going here that'll give you a somewhere near legitimate pro- 
duction rate, is one that involves low-lying transitions in the ground- rotational 
state, and that's HCN at this point. 

So in the future we can choose our molecviles particularly for production 
rate calculations, it would seem to me. In terms of comparing the interstellar 
medium to the comet, of course we get varying numbers as we go around the 
interstellar medium, but if we pick the Orion cloud, for example, and try to 
compare projected densities along a given length, we have to say that the mole- 
cules along a given length for the heart of the Orion nebula is probably about 
10^'^ integrated column density per centimeter squared. 

Now, when we try to compare that with water, we run into exactly the same 
problem we're running into in the comet, and that is that the water is masering. 
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it's time-varying, and in Orion it has about a dozen velocity components. So you 
cannot use a direct measurement of water in Orion, but you can say, "If I look 
around at the other polyatomic molecules that are seen in Orion, and I look at 
the diatomic molecules that are seen in Orion, quite typically for CO, which is 
very bright in Orion — it's about 10^ * per centimeter squared. 

And OH, from what can be deduced from absorption measurements, is about 
10*’ per centimeter squared. HCN is about 10*^, and then when we drop down 
to isocyanic acid and molecules of this complexity then we're talking about 10 * ^ , 
approximately, per centimeter squared. I'm a little hesitant on this number, 
because we just finished this detection. I haven't really completed all the anal- 
ysis on it. 

So let me raise this to 10* and again, dimethyl ether — in Orion, again is 
about 10 *^ per centimeter squared. 

Okay. So if you bootleg it like that, and say — I don't know of any peculiar- 
ity in the interstellar chemistry, given enough hydrogen and given enough oxygen, 
that would cause water to deviate very much from this scheme — then I think that 
probably within an order of magnitude either way, water should be about 10 *®, 
if you could see thermal water in Orion. You can't, so we don't have a better 
number than that, and if you wish, for your purposes you can choose 10 * ® , but 
I don't think you can choose 10* * at this point. 

B. Donn; This is a point I was going to discuss in my paper, on Wednes- 
day, on the origin of comets. It seems to me this is an important way of look- 
ing at the origin. One can, in a sense, divide the composition of comets, in- 
cluding new comets, into a number of different classes, depending on the ratios 
CO^ and Nj : CN, Cj, C 3 , NH, etc: continuum and try to fit this into some 
formation pattern which I will do on Wednesday. 



